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Abstract 

Superscaling approximation (SuSA) predictions to neutrino-induced 
charged-current charged pion production in the A-resonance region are ex- 
plored under MiniBooNE experimental conditions. The results obtained 
lO ■ within SuSA for the flux-averaged double-differential cross sections of the 

C<*j \ 7r + production for the v^ + CH2 reaction as a function of the muon kinetic 

energy and of the scattering angle, the cross sections averaged over the angle, 
the total cross section for the tt + production, as well as CCl7r + to CCQE 
cross section ratio are compared with the corresponding MiniBooNE exper- 
imental data. The SuSA predictions are in good agreement with data on 
neutrino flux average cross-sections, but a somewhat different dependence 
on the neutrino energy is predicted than the one resulting from the experi- 
mental analysis. 
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1. Introduction 

The properties of neutrinos, particularly the parameters of their oscilla- 
tions, are being studied with increasing interest as these may carry important 
information about the limits of the Standard Model. In most neutrino exper- 
iments, the interactions of the neutrinos occur with nucleons bound in nu- 
clei. The influence of nucleon-nucleon interactions on the response of nuclei 
to neutrino probes must then be considered, ideally in a model independent 
way. Model predictions for these reactions involve many different effects such 
as nuclear correlations, interactions in the final state, possible modification 
of the nucleon properties inside the nuclear medium, that presently cannot 
be computed in an unambiguous and precise way. This is particularly true 
for the channels where neutrino interactions take place by means of excita- 
tion of a nucleon resonance and ulterior production of mesons. The data on 
neutrino-induced charged-current (CC) charged and neutral pion production 
cross sections on mineral oil recently released by the MiniBooNE collabo- 
ration [1] provides an unprecedented opportunity to carry out a systematic 
study of double differential cross section of the processes, v^p — > u.~pn + , 
v^n — >■ fi~nir + , v^n — > fi~pir°, averaged over the neutrino flux. 

One way of avoiding mo del- dependencies is to use the nuclear response 
to other leptonic probes, such as electrons, under similar conditions to the 
neutrino experiments. Thus, in this paper we compare the SuSA predic- 
tions for neutrino-induced CC charged pion production cross sections with 
MiniBooNE data Jl]. The extensive analyses of scaling [2|,|3|,|4j and superscal- 
ing [5|, [g, lD, la, l9|, [l0| phenomena observed in electron-nucleus scattering lead 
to the use of the scaling function directly extracted from (e, e') data to predict 



neutrino (antineutrino)-nucleus cross sections [ll[, not relying on a particular 
nuclear structure model. Within SuSA a "superscaling function" f(ip) is built 
by factoring-out the single-nucleon content off the double-differential cross 
section and plotting the remaining nuclear response versus a scaling variable 
i/j(q,u). Approximate scaling of the first kind, i.e., no explicit dependence of 
f(ip) on the momentum transfer q, can be seen at transfer energies below the 
quasielastic (QE) peak. Scaling of second kind, i.e., no dependence of f(ip) 
on the mass number, turns out to be excellent in the same region. When 
scaling of both first and second types occur, one says that superscaling takes 
place. 

The analyses of the world data on inclusive electron- nucleus scattering [7| 
confirmed the observation of superscaling and thus justified the extraction 



of a universal nuclear response to be also used for weak interacting probes. 
However, while there is a number of theoretical models that exhibit superscal- 
ing, such as for instance the relativistic Fermi gas (RFG) p, |6j, the nuclear 
response departs from the one derived from the experimental data. This 
showed the necessity to consider more complex dynamical pictures of finite 
nuclear systems - beyond the RFG - in order to describe the nuclear response 
at intermediate energies. SuSA predictions are based on the phenomeno- 
logical superscaling function extracted from the world data on quasielastic 



electron scattering [12|. The model has been applied to neutral current scat 



tering [13[ and it has also been extended to the A-resonance region [TJJ where 
the response of the nuclear system proceeds through excitation of internal 
nucleonic degrees of freedom. Indeed, a non-quasielastic cross section for 
the excitation region in which nucleon excitations, particularly the A, play 
a major role was obtained by subtracting from the data QE-equivalent cross 
sections given by SuSA [9|, ll4| . This procedure has been possible due to the 
large amount of available high-quality data of inelastic electron scattering 
cross sections on 12 C, including also separate information on the longitudi- 
nal and transverse responses, the latter containing important contributions 
introduced by effects beyond the impulse approximation (non- nucleonic). 
The SuSA procedure has been already employed to describe the non pio- 



nic (QE) cross-section of the MiniBooNE v- and //-nucleus cross-section [15|, 



16l . Il7| . Here we extend the analysis to CC pion production cross-section 
measured at MiniBooNe, that from the theoretical point of view can be seen 
as more challenging. For instance, A properties in the nuclear medium, as 
well as both coherent and incoherent pion production for the nucleus should 
be considered in any theoretical approach, while in the SuSA procedure they 
are included phenomenologically extracted from the electron scattering data. 
All what is assumed within SuSA approach is an internal factorization of the 
nuclear response to a weakly interacting probe into a single-nucleon part 
and a 'nuclear function' accounting for the overall interaction among nucle- 
ons. As mentioned before, the SuSA assumptions have been tested against 
a great deal of electron-nucleus scattering data with fair success. The fac- 
torization assumption allows to apply the same nuclear responses derived 
from electron scattering to neutrino-induced reactions, with a mere use of 
the adequate single-nucleon terms for this case. To show the importance 
of nuclear interaction effects as predicted within SuSA, as a reference, we 
also show results obtained within the RFG, with no interactions among nu- 
cleons, for which the scaling function in the A-domain is simply given as 



3 

/^p G (^a) = -7(1— ''Pa 2 ) 9(1 — ip a 2 ) with ip& the dimensionless scaling variable 

extracted from the RFG analysis that incorporates the typical momentum 
scale for the selected nucleus |S|, |ll . 



In Fig. [T] we compare the A-region SuSA [UJ] and RFG scaling functions, 
that we use in our study. Here the data refer to 12 C and 16 and span a 
large range of energies (from 0.3 to 4 GeV) and scattering angles (from 12 to 
145 degrees). The experimental points in Fig. [1] are extracted by subtracting 
from the total cross sections the quasielastic contribution calculated using 
the universal QE scaling function /qe(V>)- This analysis does not include 
the contribution associated to meson-exchange currents (MEC), which are 
important in the region between the QE and A peaks and are responsible of 
the disagreement between the data and the fit (red curve) at large negative 
values of ipA- These currents are mediated by virtual pions and do not corre- 
spond to the production of real pions; hence they should be included in the 
"quasielastic" response, which for the MiniBooNE experiment corresponds to 
the absence of real pions in the final state but not to the pionic data discussed 
in this work. The contribution of the MEC to the QE neutrino and antineu- 
trino scattering has been evaluated in the SuSA framework in two recent 



papers [15|, [17| . We also note that the fit shown in Fig. [TJ is restricted only to 
excitation energies at and below the A resonance peak, where the response is 
dominated by the A; at higher energies other resonances and eventually the 
tail of deep inelastic scattering (DIS) contribute. This explains the difference 
between the phenomenological superscaling function, which aims to describe 
the A resonance peak, and the data observed in Fig. CD at positive ^a- values. 

2. Formalism 

In what follows we present the results of applying the SuSA and RFG A- 
scaling function to neutrino-induced CC charged pion production. We follow 



the formalism given in [ll]. The charged current neutrino cross section in 



the target laboratory frame is given in the form 



d 2 a (Gcos6 c k') 2 ( |Q 2 , 



dD.dk! 2?r 2 V 4 ee' 

where Q, k' and e' are the scattering angle, momentum and energy of the 
outgoing muon, G is the Fermi constant and 9 C is the Cabibbo angle. The 
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Figure 1: (Color online) The SuSA scaling function in the A-region / a (V>a) (solid line) 
extracted from the world data on electron scattering ll|. The dotted line shows the 
scaling functions / a (V'a) in the RFG model. 



function J 72 depends on the nuclear structure through the R responses and 
can be written as 11, 1 



J 7 — Vcc-Rcc + 2Vcl-Rcl + Vll-Rll + Vt-Rt + 2V T 'R T i 

that is, as a generalized Rosenbluth decomposition having charge-charge 
(CC), charge- longitudinal (CL), longitudinal-longitudinal (LL) and two types 
of transverse (T,T') responses (-R's) with the corresponding leptonic kinemat- 
ical factors (V's). The nuclear response functions in A-region are expressed 
in terms of the nuclear tensor W^ u in the corresponding region. The basic ex- 
pressions used to calculate the single-nucleon cross sections are given in 
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These involve the leptonic and hadronic tensors as well as the response and 
structure functions for single nucleons. A convenient parametrization of the 
single-nucleon W + n — > A + vertex is given in terms of eight form-factors: four 
vector (C34 5 6 ) and four axial (C^ 4 5 6 ) ones. Vector form factors have been 
determined from the analysis of photo and electro-production data, mostly 
on a deuteron target. Among the axial form factors, the most important 
contribution comes from C§. The factor C$, whose contribution to the dif- 
ferential cross section vanishes for massless leptons, can be related to C^ by 
PCAC. Since there are no other theoretical constraints for C^ 45 (g 2 ), they 
have to be fitted to data. We use two different parameterizations: the one 



given in 19] where deuteron effects were evaluated (authors estimated that 



the latter reduce the cross section by 10%), denoted as "PR1", and the one 



from [20[, called "PR2". 

With these ingredients, we evaluate the cross section for CC A ++ and 
A + production on proton and neutron, respectively. Once produced, the A 
decays into irN pairs. For the amplitudes A of pion production the following 
isospin decomposition applies: AiyiV ~~ >* l~pK + ) = A3, Aiyin — > l~nn + ) = 
lA 3 +^p- A±, A(uin ->■ l~pTi Q ) = -^As+^At, with .A3 being the amplitude 
for the isospin 3/2 state of the ttN system, predominantly A, and A\ the 
amplitude for the isospin 1/2 state that is not considered here. 

3. Results and discussion 

First we present RFG and SuSA predictions for the double-differential 
cross section for CC neutrino-induced n + production on CH 2 averaged over 
the neutrino flux $(e^), namely 



d 2 a 



dT^d cos 9 & t 



d 2 a 



dT^d cos 9 



®{e„)de u , (2) 



where T M and 9 are correspondingly the kinetic energy and scattering angle 
of the outgoing muon, e u is the neutrino energy and $ tot is the total inte- 
grated v^ flux factor for the MiniBooNE experiment ($tot = 5.19 x 10~ 10 
[z/^/cm 2 /POT]). The double-differential cross section averaged over the neu- 
trino energy flux as a function of the muon kinetic energy T M is presented in 
Fig. [2j Each panel corresponds to a bin of cos 9. The PR2 parametrization 
has been considered. Results with the PR1 parameterization are about 5% 
higher, that is a measure of the degree of uncertainty that we expect from 
the choice of the single-nucleon response for this reaction. We compare the 
predictions of SuSA and RFG with the MiniBooNE data [l|]. The nuclear 
target has been considered as carbon and hydrogen in the mineral oil target. 
Fig. [3] shows SuSA and RFG predictions for the double-differential cross sec- 
tion averaged over the neutrino energv flux as a function of the scattering 
angle at fixed T M compared with data [l[. 

Figs. |2] and [3] show a good agreement between data and the SuSA pre- 
dictions for the flux-averaged double-differential cross sections. This applies 
to both parameterizations of the vector and axial form factors. RFG results 
have similar shape as SuSA ones but, as expected, they overestimate the 
data to a large extent. 
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In Fig. H] are shown the results obtained by integrating the flux-averaged 
double-differential cross sections over angle: 

w^^s^iys^)^^ < 3) 

The total cross section for ir + production as a function of the neutrino 
energy along with the MiniBooNE data are displayed in Fig. [5j Poorer 
agreement with data than for the flux-averaged cross sections presented in 
Figs. [2HH is clearly observed. The data seem to follow a more linear depen- 
dence with the energy up to 2 GeV than the theory. However, before drawing 
definite conclusions, one has to consider on one side that the nuclear response 
extracted from the phenomenological electron data includes the whole (real) 
pion production strength (virtual pion contribution via MEC has already 
been removed) while in MiniBoone the data are sensitive only to the cases 
where a real pion is seen away from the nucleus. This means that if the real 
pion is produced and then absorbed during the final state, it will not be seen 
in the observed MiniBoone pionic cross-sections. On the other hand, the un- 
folding procedure used to extract the data of Fig. [5] is to some extent model 
dependent. Thus these data are less direct and we consider the comparison 
with the data of Figs. [2H3] to be of more significance. It is worth mentioning 
some recent publications where the problems with the reconstruction of neu- 



trino energy E v are discussed in more details [211. I22L [23| and also the review 
by H. Gallagher, G. Garvey, and G.P. Zeller [24|, where the authors con- 
sider in depth neutrino-nucleus interactions (in the medium-energy regime, 
0(1 GeV)) in respect to the modern neutrino oscillation experiments. 

Figure |6] shows the ratio of CCl7r + (CC single-pion production) to CCQE 
(CC quasielastic scattering) cross sections from SuSA, SuSA+MEC (2p-2h 
meson-exchange current) [l6|j , and RFG approaches in comparison with the 
MiniBooNE data corrected for final state interactions. All these ratios have 
been rescaled to an isoscalar target [25|. The results are obtained on the 



basis of total cross sections for CCl7r + (given in Fig. [5]) and CCQE 16(. A 
similar conclusion as the one in the previous figure could be drawn here. It 
seems that there is too much n production strength below 1.2 GeV, and too 
little beyond that, compared to data. 

Before concluding we would like to remind that the A scaling function 
used in our approach represents a good fit of electron scattering data only 
in the region below the A-peak, as shown in Fig. [TJ As anticipated, higher 



resonances and the tail of DIS come into play at high momentum and energy 
transfer. However, in the kinematical conditions of the MiniBooNE experi- 
ment the dominant contribution is associated to A-excitation and therefore 
we only include the latter in our analysis. This is supported by the fact 
that inclusive electron scattering data in the same kinematical domain are 
reasonably well reproduced by using the pure A-superscaling function (see 
Ref. H). " 

4. Conclusion 

Summarizing, in this paper we present results for the cross sections 
of neutrino-induced CC n + production obtained with the SuSA and RFG 
(shown as reference) models. The SuSA approach provides neutrino-nucleus 
cross-section predictions, based on the observed nuclear response to elec- 
tron projectile and the universal character of the scaling function. Notice 
that SuSA predictions incorporate effects of final state interaction (FSI), the 
properties of the A resonance in the nuclear medium, both the contribu- 
tion of coherent and incoherent production, etc.. The role of the FSI on 
the one-pion production has been considered for instance within the GIBUU 



transport model [26|, where it was shown that in order to reproduce the data, 
the total cross section obtained with FSI included has to be multiplied by 
a factor of 1.5. Here we show that SuSA predictions are in good agreement 
with the MiniBooNE experimental data for pionic cross-section in the case 
of the flux averaged data, while some disagreement remains in the compar- 
ison to unfolded neutrino energy data. Notice that the accordance between 
SuSA and data here is better than the one for the non-pionic case, where the 
model was found to underpredict the data unless meson exchange currents 
were explicitly included [16[. We conclude that the SuSA scaling function 
for the A-region (extracted from electron scattering experiments) and its ex- 
tension to neutrino processes may be very useful to predict cross sections for 
neutrino-induced CC tt + production, not relying on specific models. 
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Figure 2: (Color online) The double-differential cross section averaged over the neutrino 
energy flux as a function of the muon kinetic energy T^ obtained by SuSA (solid line) 
and RFG (dotted line) A-rcgion scaling functions. In each subfigure the results have been 
averaged over the corresponding angular bin of cos 6. "PR2" parametrization [20( is used. 
The results are compared with the MiniBooNE data [l|. 
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Figure 3: (Color online) Same as Fig. [5] but the double-differential cross section averaged 
over the neutrino energy flux as a function of the scattering angle, are presented. 
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Figure 5: (Color online) The total cross section for tt + production are compared with the 
MiniBooNE data [l|. For vector and axial form- factors two parameterizations, "PR1" [19( 
and "PR2" [H, are used. 
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Figure 6: (Color online) The results for CCi7r + to CCQE cross section ratio are compared 
with MiniBooNE data (corrected for final state interactions and rescaled for an isoscalar 
target) J25|. 
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